Journal of Ecology

Coflowering, the temporal overlap in flowering among plant species, can influence plant fitness through its effects on heterospecific pollen transfer and competition for pollinators and resources. This temporal overlap is changing as plant species respond differently to climate change. However, difficulties in quantifying and comparing patterns of coflowering across diverse, multispecies communities have limited progress in understanding climate-driven shifts in coflowering. Network-based approaches offer a promising solution to this limitation. Here, we investigate how warming and altered precipitation influence coflowering network structure, flowering phenology, and seed production in 12 annual California serpentine grassland plant species in a mesocosm experiment, using coflowering network analysis and structural equation modeling. Communities consistently grouped into two phenological modules corresponding to early- and late-season annuals. Warming reduced the overall amount of coflowering in the community, driven primarily by weakened and less diverse coflowering among late-season species. Soil moisture moderated these effects: late-season coflowering was strongest under cool-wet conditions and weakest under warm conditions regardless of soil moisture. In contrast, early-season species had a more fixed phenological window, and maintained stable coflowering relationships across climates. Growth form (grass or forb) and origin (native vs. non-native) did not predict coflowering responses. We additionally found that temperature and soil moisture influenced seed production largely through their effects on flowering phenology, with early- and late-season species exhibiting distinct phenological responses and sensitivities. Our results demonstrate that warming and precipitation change can spread out flowering times within functional groups, reshaping community-wide coflowering networks through species-specific responses to altered climate conditions.

Repeated grazing and trampling can reduce forb populations in South African mesic grasslands by limiting photosynthesis and depleting carbohydrate reserves stored in underground organs. We examined the effects of repeated defoliation on growth, bulb starch reserves and water status in the geophyte Eucomis autumnalis over 468 days. Four intense summer and autumn defoliations reduced above-ground production, bulb mass, starch reserves and bulb water pools, with the strongest impacts occurring in late summer and autumn when plants normally replenish reserves. Repeated defoliation disrupted seasonal source sink dynamics, forcing resources into regrowth and driving bulbs towards energetic bankruptcy. These results highlight the importance of avoiding heavy autumn grazing and providing periodic full-year rests to conserve geophytic forbs in mesic grasslands.

AO_SCPLOWBSTRACTC_SCPLOWMountain gradients facilitate our understanding of species range limits, competition dynamics, stress-resilience trade-offs, and determinants of vegetation zone boundaries. Forest compositional models often use altitude as the main predictor, a proxy for temperature that is defensible where floristic transitions are gradual and climate relationships are linear. However, mountains with distinct assemblages, representing tropical gradients or areas with complex biogeographic history, require a modeling framework that reflects non-linear dynamics or interactions between environmental factors, including outlier events (rather than mean conditions). Our study system encompasses both tropical and temperate forests along a broad ([~]3000 m) altitudinal gradient, positioned within a narrow latitudinal band (< 1{degrees}) and composed of mature, continuous forest in the Bhutan Himalaya. To represent the breadth of climatic conditions experienced over a trees lifetime, we used a Bayesian modeling paradigm and integrated multi-generational field knowledge to develop a priori hypotheses and informed priors, with consideration of monsoon seasonality and possible ecophysiological thresholds. Our approach followed three stages (the Pattern, the Mechanism, the Test). Specifically, we interpolated microclimate data and derived custom metrics based on thermodynamics, propagating uncertainty into subsequent models to test whether climate posteriors outperformed altitude in explaining growth form partitioning. For spatial patterns, we identified six distinct vegetation zones (encompassing 145 species from 57 families), with a mid-gradient peak in richness at the tropical-temperate transition zone, and convergence of deciduousness at either end of the gradient. For individual growth forms, abundance was tied to different ecological mechanisms, explained by adaptations to climatic stressors and competition trade-offs. For instance, evergreen broad-leaved dominance was linked to ephemeral cloud immersion, whereas tropical deciduous species were affiliated with higher vapor pressure deficit at lower altitudes. Most importantly, compositional (between-group) models showed that the interaction between frost events and fog probability (air saturation prior to the dry season) governed growth form partitioning more than any single factor; temperate deciduous species, confined to a narrow altitudinal band, exemplified this finding. Our methodological approach is transferable to other data-sparse mountain systems, and our results highlight the vulnerability of unique habitat types and montane endemics under climate change scenarios that alter the fog-frost dynamics. Second abstract in DzongkhaTo see the second abstract in Dzongkha, the official language of Bhutan, please visit our Zenodo site: https://doi.org/10.5281/zenodo.19081441.

A central question in biodiversity conservation is whether species can sustain viable population under current and future atmospheric N deposition. Assessing species viability under N deposition requires demographic studies integrating species vital rates responses to long-term N deposition across different levels. However, studies of this nature are rare. Our integral projection models (IPMs), parameterized with demographic data, revealed differing responses of two functionally similar coexisting species, Stipa bungeana and Leymus secalinus, to 12 years of N deposition at low N addition levels (1.15 and 2.30 g N m-2 yr-1) and high N addition levels (4.60, 9.20, and 13.80 g N m-2 yr-1) on the Loess Plateau grasslands. We found that the reduced survival across N addition levels was partially compensated by increased contributions from growth, shrinkage, and fecundity, alleviating the population decline of S. bungeana (with a longer lifespan and generation time) under different N additions. Contrasting, more positive correlations among vital rate enabled the population of L. secalinus (with a shorter lifespan and generation time) to track N additions, with population growth under low N additions and population decline under high N additions. Our results illustrate that the demographic response to N deposition may vary considerably between functionally similar coexisting species, and species with demographic compensation can buffer populations against N deposition while with demographic lability enable populations to track N deposition. Furthermore, our study demonstrates the potential of using life-history traits to predict species viability under N deposition, thereby informing biodiversity conservation under global change.

While mutualistic symbioses with nitrogen-fixing bacteria enable plants to access fixed nitrogen, they also require substantial carbon investment. Under carbon limitation, such as shading, shifts in biomass allocation can decouple symbiotic investment from leaf and root growth, potentially compromising plant nitrogen status. Because shading shifts biomass allocation toward light acquisition, it could influence nitrogen fixing symbiosis in two opposing ways. If nodulation remains coupled to leaves rather than roots, nitrogen status should be maintained despite reduced root growth. Alternatively, if root growth constrains nodulation, nitrogen status should decline. We tested these hypotheses by manipulating light availability (full sunlight vs. 50% shade) and quantifying biomass allocation and symbiotic nodulation. Under shading, plants allocated proportionally more biomass to shoots than to roots and invested less biomass in root nodules. Relationships between nodulation and leaf or root biomass differed between treatments but converged with increasing plant size, although shaded plants never attained the root biomass observed in full sunlight. Leaf nitrogen concentration was maintained under shading because nodulation remained coupled to leaf investment despite reduced root allocation. These findings highlight that, under carbon limitation, maintaining leaf and nodule coupling enables plants to reduce nodule investment without compromising the nitrogen benefits of symbiosis.

The increasing frequency and intensity of heatwaves under climate change highlight the need to understand how biodiversity regulates forest canopy thermal dynamics. Although tree diversity can buffer the microclimate, its effects on canopy temperature and the role of mycorrhizal symbioses remain unclear. We addressed this question in the MyDiv tree diversity experiment in Germany, where tree species richness (1-, 2-, and 4-species mixtures) and mycorrhizal types (arbuscular, ectomycorrhizal, and mixed) are factorially manipulated. During the 2024 growing season, we conducted nine uncrewed aerial vehicle (UAV) surveys using integrated thermal and LiDAR sensors to quantify canopy temperature and structural complexity, together with measurements of leaf water content, specific leaf area, soil moisture, and vapour pressure deficit (VPD). Increasing tree diversity generally reduced canopy temperature, although the strength of this relationship varied seasonally and among mycorrhizal types. Cooling effects were strongest during peak summer heat and were more pronounced in arbuscular mycorrhizal (AM) communities than in ectomycorrhizal (EM) and mixed (AM+EM) communities. In contrast, EM communities exhibited greater canopy structural complexity, whereas AM communities maintained higher soil and leaf water content. Structural complexity increased with tree diversity but did not necessarily result in greater canopy cooling. Structural equation modelling revealed that forest thermal buffering emerged through complementary structural and hydraulic pathways, whose relative importance shifted seasonally, with hydraulic regulation becoming increasingly important under hotter and drier conditions. By linking canopy temperature, canopy structure, and plant water relations, our study provides mechanistic insights for understanding how multiple facets of biodiversity regulate forest thermal buffering under climate warming.

Artificial light at night (ALAN) can alter the movement of nocturnal pollinators such as moths, with potential consequences for pollination services. However, most studies have focused on the presence or properties of lights, while the role of spatial lighting configuration remains poorly understood. We conducted a small-scale field experiment to test how different lighting configurations affect pollination success in moth-pollinated plants, using low-intensity LED garden lights. Potted phytometer plants of three moth-pollinated species were exposed to one of three treatments: no-light control, isolated LED point lights, or 25 m linear arrays of multiple LED lights. We quantified pollination success as both the probability of seed set and reproductive output through seed and capsule production. Experimental lighting affected plant reproductive success, with light arrays increasing seed-set probability, seed and seed capsule number relative to unlit controls, whereas point lights showed no clear effect. These results provide preliminary evidence that the spatial configuration of artificial lights may influence nocturnal pollination outcomes. Future work combining phytometer assays with direct tracking of moth movement is needed to assess whether light arrays facilitate or redirect the dispersal of nocturnal pollinators.

O_LISpecies interactions are increasingly recognised as temporally dynamic. For plant-pollinator networks, evidence shows that interactions vary not only seasonally but also over the diel cycle. However, we still know little about what determines their diel structure, whether this structure is important for robustness to species loss, and which plant traits are associated with the roles plants play in this structure. These gaps are fundamental, because knowing what shapes networks over the diel cycle is required to predict the effects of global change drivers, such as light pollution, that may shift the timing of interactions. C_LIO_LIUsing 22 plant-pollinator networks sampled over morning, afternoon, and night, we addressed these gaps by applying a multilayer framework to characterise their diel structure, link it to robustness, and test which plant traits are associated with plant roles across diel periods (participation, versatility) and within them (centrality). C_LIO_LIDiel structure was non-random: interactions were segregated among diel periods yet integrated through plants visited across the diel cycle. Networks were more robust to simulated species loss when interactions were on average more evenly distributed across periods and when plants were more strongly interconnected among periods, although this benefit diminished when both properties were high simultaneously. The association between plant traits and their roles shifted with the temporal scale. Across diel periods, structural traits were the stronger predictors: taller plants were visited more evenly across the diel cycle (higher participation), whereas shorter plants shared pollinators with plants from multiple periods (higher versatility), potentially mediating indirect effects among them. Within diel periods, floral visual cues became more influential, with achromatic contrast the most consistent predictor: at night, plants with brighter flowers were well visited within the period, sharing pollinators with other plants of that period (higher centrality). C_LIO_LIThese findings identify diel structure as a functional axis of network organisation and indicate that plant-pollinator networks are assembled hierarchically: structural traits set a baseline across the diel cycle, whereas light conditions determine which traits matter within periods, ultimately defining distinct temporal pathways vulnerable to global change. C_LI

SummaryO_LIIn woody perennials, reproductive allocation considered as a fraction of NPP (RA) has been rarely quantified, especially tracked across plant lifetimes, measured in biologically meaningful currencies, and separated from standing biomass. C_LIO_LIWe addressed this gap by measuring dry mass, nitrogen, and phosphorus for every aboveground tissue type for 14 iteroparous perennial shrubs tracked across their full lifetimes, calculating RA using four accounting schemes and three currencies. C_LIO_LIRA was substantially higher by mid-life than typical estimates from ecosystem-scale studies. Distinguishing standing biomass from yearly production further revised RA upwards. Using a nutrient currency (especially P) increased the perceived costs of reproduction. Propagules were highly nutrient-enriched and reproductive accessory costs consumed more nutrients than did the propagules themselves. Considering N and P resorption from senescing leaves and wood shrank the effective vegetative nutrient budget, further concentrating net annual nutrient demand in reproductive tissues. C_LIO_LIOur results highlight high investment in RA for woody perennials, especially using nutrient currencies. Broadly similar allocation patterns were observed across species with different functional traits and lifespans, suggesting generality that may apply across biomes. Widespread underestimation of RA in forest growth models likely overestimates the proportion of NPP available for vegetative growth, leading to substantial errors in predictions. C_LI

A central question in biodiversity conservation under global change is whether species maintain viable populations under both mammal disturbance and climate warming. This requires demographic studies that integrating vital rates responses to mammal disturbance and climate warming across an entire life to. Using Integral Projection Models parameterized with demographic data, we found the population growth rates of Thermopsis lanceolata under both ambient and warming conditions, initially decreased on new mounds, further declined on seminew mounds, but eventually exceeded initial levels on old mounds. This stage-dependent responses were largely driven by clonal reproduction (i.e., clonal production and/or ramet size distribution), which emerged as both the most sensitive vital rate and the primary contributor to variation in population growth rates across recovery stages. Additionally, we found that the combined effects of plateau zokor disturbance and warming on population growth rates of new mounds was greater than the sum of their individual effects, leading to population decline on new mounds. Such synergistical effects was mainly due to a larger decrease of ramet size distribution. These findings suggest that multifactorial experiments are important for biodiversity research, rather than merely adding single effects on population dynamics. In addition, clonal reproduction may be a key vital rate for population maintenance under global change.

Predicting species interactions remains a major challenge, as multiple species attributes operate simultaneously and their relative importance may vary seasonally and with temporal resolution. Here, we assess how the relative contributions of abundance, trait matching, and phenology to plant-pollinator interactions vary through a flowering season and across temporal resolutions using interaction data from three European botanical gardens. We show that predictive models explain a substantial proportion of variation in visitation patterns, with floral and pollinator abundances consistently explaining most variation across the flowering season. Trait matching between pollinator body size and floral size also contributes to visitation patterns, playing a secondary but persistent role in shaping interactions, while phenology plays a relatively minor role at broader temporal resolutions but becomes more important at finer temporal scales. Additionally, null models accounting for spatio-temporal variation in floral and pollinator abundance reveal consistent patterns of pollinator preference and avoidance, indicating that abundance alone cannot explain the observed interaction patterns. Our results show that seasonal variation and temporal resolution differentially shape the importance of species attributes, highlighting the need for multi-variable approaches that account for temporal dynamics to accurately explain and predict ecological interactions.

Climate change is accelerating species losses in ecosystems across the world. Seed germination is a critical, climate-dependent phase of the plant life cycle; however, the ecological determinants of germination climate niches within diverse landscapes and across functional types (FTs) are still not well understood. In this study, we characterized seed germination temperature and water availability niches for 28 species that represent different FTs (tree, shrub, grass, forb) and vegetation types (grassy woodland, dry and wet forests) within a temperate bioregion (Sydney, Australia). We tested whether ecological determinants, specifically species climate of origin, seed traits, FT and vegetation type explain germination niches and predicted spatial and temporal patterns of germination potential across the landscape under high and low emission scenarios. We found wide variation in thermal and hydric germination niches among species. Optimal germination temperature (thermal niche) was predicted by FT, climate of origin and seed traits, such that shrubs, cool-origin species, and species with large seeds had significantly cooler optimal temperatures for germination. We also quantified spatial and temporal changes in germination potential to identify vulnerable areas and FTs. We found strong species-specific seasonal patterns in germination potential with future climate shifts affecting FTs differently; germination of woody species declined more than forbs. Future germination potential was predicted by historical climatic conditions, with warmer and drier localities being more vulnerable. Overall, our findings demonstrate that species germination responses to climate change depend on FT, seed traits, and species climate of origin, with woody species and warmer, drier parts of the landscape emerging as being particularly vulnerable to declines in recruitment. Our study provides a mechanistic understanding of germination responses to temperature and water availability, enabling predictions of vulnerable species and areas for conservation under climate change, and inform large-scale ecosystem restoration approaches through improved species selection and sowing times.

O_LIForests are vertically structured ecosystems where light attenuation, microclimatic buffering and resource availability occur along continuous gradients from the canopy to the understorey layers. Despite this complexity, studies rarely integrate trophic interactions across vertical layers, overlooking how vertical forest structure shapes consumer communities through abiotic and biotic pathways. C_LIO_LIIn this study, we combined layer-specific measurements of plant diversity, structure and microclimate with arthropod sampling in the canopy and understorey, as well as bird survey data, in a temperate forest. Applying Bayesian structural equation models with explicitly defined directional pathways, we modelled both consumer biomass and abundance across trophic levels and vertical layers. C_LIO_LIAbundance measures were predominantly filtered by local layer conditions, while biomass responded to conditions across layers, reflecting stand-level energy flow. This suggests that these two metrics capture fundamentally different ecological processes. Canopy conditions consistently predicted understorey arthropod abundance across trophic levels, while the reverse was not observed, suggesting a strong asymmetric downward propagation of canopy-driven effects. Furthermore, trophic interactions between arthropod primary and secondary consumers remained largely stratified within vertical layers, suggesting a vertical food web compartmentalisation rarely shown in structurally complex aboveground systems. C_LIO_LIPlant diversity, structure and microclimate shaped consumer communities mainly through the modulation of resource availability and plant apparency, with effects varying across vertical layers, trophic levels and taxonomic groups. Through complementary mechanisms, plant diversity likely determined the variability of resources available to consumers at different trophic levels. Structural properties, in contrast, potentially drove the spatial redistribution of these resources through light attenuation and microclimatic buffering, which in turn influenced the physiological capacity of consumers to access and exploit available resources. C_LIO_LIThese findings demonstrate that vertical stratification mediates trophic pathways in a highly directional manner, with canopy characteristics playing a disproportionate role in structuring the forest community across layers. Integrating layer-specific structural and trophic indicators into forest biodiversity assessments and management strategies is therefore essential to fully evaluate biodiversity dynamics and multifunctionality in structurally complex forest ecosystems. C_LI

Grasslands are vulnerable to increasing drought with global warming, but process-based models lack the mechanistic knowledge required to predict the magnitude of drought impacts. While a plant hydraulics framework has been successful in advancing process understanding of drought responses in trees, and how drought responses vary across rainfall gradients, similar approaches have rarely been applied to grasses. Here, we quantified the progression of key drought response processes in sixteen dominant perennial grasses (seven C3 and nine C4) with differing climatic origins across eastern Australia. We found that stomatal closure, hydraulic impairment and leaf browning occurred concurrently, in contrast to the progressive sequence typically observed in trees. We also found that drought response traits were not correlated with species climate of origin. The early impairment of leaf hydraulic conductance and leaf browning along with the lack of correlation with climate of origin suggest that grasses may employ fundamentally different strategies to adapt to low water availability than trees. These results highlight the need for grass-specific parameterization of drought responses in process-based models.

A central question in restoring degraded grasslands is whether grazing removal can sustain viable plant populations under both current and future warming conditions. Addressing this question requires demographic studies integrating vital rates responses to grazing removal and climate warming throughout a species life cycle. However, studies of this nature are rare. Using stochastic Integral Projection Models parameterized with four years (2020-2023) of demographic data, we find that nine years of grazing removal increases the stochastic population growth rate (log{lambda}S) of two coexisting herbaceous plants, Carex atrofusca and Sibirotrisetum sibiricum at two altitudes (3,700 m and 4,000 m) in an alpine grassland on the Tibetan Plateau. Although individual survival declines following grazing removal, these negative effects are overcompensated by enhanced plant growth, ultimately promoting log{lambda}S in both species. However, the benefits of grazing removal are cancelled under nine years of in situ active warming (+2{square}), where no demographic compensation occurred (i.e., vital rates change in the opposite directions among populations), and log{lambda}S are even lower than those under grazing. Our findings suggest that while grazing removal is a sustainable management strategy under current climate conditions, it may not remain effective under projected warming, providing valuable information for sustainable population management under global change.

Tropical rainforests generate and maintain their own microclimate regimes and the resultant cooler, more humid and stable environments foster the hyperdiversity typical of these ecosystems. Although the temporal and spatial (horizontal) distributions of microclimates have been relatively well studied, the vertical dimension has received less attention, and little is known about how forest disturbance affects the vertical stratification of microclimates in tropical forests. In this study, we examine how the vertical distribution of temperatures varies between undisturbed and burned Amazonian forests. We installed five vertical transects with temperature dataloggers distributed at 7 different heights to collect data over multiple days during the end of the dry season. We investigated how anthropogenic disturbance (fire) mediates the vertical stratification of microclimate and whether microclimate buffering (i.e, the difference between understorey and canopy temperatures) varies according to the forest structure. We showed that anthropogenic disturbance can cause an inversion in the vertical stratification of microclimates, with burned forests having hotter temperatures (up to 2 {degrees}C) in the understorey than in the canopy during the day - the opposite of what is found in undisturbed forests (typically 3 {degrees}C cooler). During the night, while understorey and canopy temperatures are similar in undisturbed forests, we found that, in burned forests, understorey temperatures were up to 2 {degrees}C cooler than in the canopy. Microclimate buffering by day was best explained by aboveground carbon stocks, with higher temperature buffering in more carbon rich forests. Our study shows that anthropogenic disturbance alters the vertical stratification of temperatures in Amazonian forests, leading to significant temporal changes along the diel cycle. Future research should focus on understanding these changes across a wider range of disturbance regimes, and explore the consequences for biodiversity and ecosystem functions from the canopy to the forest floor.

Elevational gradients provide a framework for understanding how environmental filtering reorganises communities and interactions, but plant-pollinator interactions along temperate forest elevational gradients remain overlooked. We studied early-spring understorey communities at four forest sites spanning the foothills towards the timberline (450-1,000 m a.s.l.) in the Krkono[s]e Mountains, Czechia. Across six transects per elevation, we quantified flowering plant species richness, floral resources and traits, and video-recorded flowers, yielding 4,003 pollinator visits. We analysed elevational patterns in species richness, community composition, floral traits, and quantitative network characteristics. Visitation frequency and flowering plant and pollinator species richness peaked at intermediate elevations. The contribution of dipteran relative to hymenopteran pollinators increased towards higher elevations, principally because of non-hoverfly flies, whereas individual bee groups showed no uniform response. Floral resources and traits showed no uniform elevational responses, although total nectar sugar availability peaked at the highest site because of the dominant Vaccinium myrtillus. Most notably, both network-level specialisation and mean species-level specialisation were generally greater at the two higher elevations, whereas nestedness was lower and other network characteristics showed no consistent patterns. These findings suggest that shifts in pollinator composition and dominant floral resources potentially shaped interactions along the gradient. The increasing specialisation with elevation contrasts with the generalisation often expected under reduced partner availability and indicates that forest networks may follow elevational patterns not predicted from open habitats. Despite limited site-level replication, this study provides, to our knowledge, the first community-wide characterisation of plant-pollinator interactions along a temperate forest elevational gradient and identifies patterns requiring evaluation across replicated gradients.

Rapid urbanization is altering the seasonal life cycles, or phenology, of urban trees. Diverse species composition and heterogeneous urban thermal environments can generate pronounced phenological variation within cities, yet such fine-scale seasonality remains poorly understood. This limits the ability to anticipate localized ecosystem services and disservices, including canopy cooling and pollen exposure. Here, we analyzed individual-tree phenology and their response to temperature across 76 cities in the contiguous United States by integrating municipal street-tree inventories, PlanetScope satellite time series, and 1-km near-surface urban air-temperature data. We quantified intra-city and cross-city phenological variation, estimated species-level associations between phenology and fine-scale urban temperature, and tested whether these associations vary with regional climatic context. Intra-city phenological variation rivaled cross-city variation for 33.3% of species in spring and 53.7% in fall. Temperature-phenology associations were widespread: warmer local urban environments were associated with earlier spring onset in 51.5% of cases and delayed fall senescence in 66.5%, with approximately two-thirds of these responses cascading to longer growing seasons. These associations varied with regional climate, generally showing stronger spring advancement and fall delay in cooler and wetter cities, but weaker or reversed responses in hotter or drier cities. Together, these results suggest that urban tree phenology is shaped by interactions among fine-scale urban temperature, species identity, and regional climate. Accounting for this heterogeneity can improve predictions of urban forest function and inform climate-resilient species selection and management. Plain Language SummaryStreet trees within the same city do not all leaf out in spring or lose their leaves in fall at the same time. This seasonal timing, called phenology, affects services such as shade and cooling, as well as disservices such as pollen exposure. We studied individual tree phenology and their temperature response in 76 U.S. cities using street-tree records, satellite observations, and local air-temperature data. We asked how much phenology varies within cities, how it responds to neighborhood-scale temperature differences, and whether these responses depend on regional climate. Variation within a city was often as large as variation among cities, especially in fall. Trees in warmer parts of cities generally leafed out earlier and lost their leaves later, often extending the growing season. However, these responses were stronger in cooler and wetter cities and weaker or sometimes reversed in hotter or drier cities. These findings show that urban tree seasonality depends on local temperature, species identity, and regional climate, which should be considered when selecting and managing trees for future cities. Key PointsO_LISubstantial within-city phenological variation is not captured by city-level averages. C_LIO_LIFine-scale temperature shaped phenology and growing-season length, but responses varied across species, cities, and phases. C_LIO_LITemperature sensitivity depended on regional climate, weakening or reversing in hotter or drier cities. C_LI

The climatic transition from Marine Isotope Stage 3 (MIS3) to the Last Glacial Maximum (LGM) had caused widespread vegetation change. Despite the dynamic equilibrium between vegetation and climate, the specific role of functional composition in vegetation response to climate change was inadequately understood. Here, we analyzed the long-term trajectories of palynological diversity, vegetation coverage and community-weighted-mean (CWM) functional traits based on EH22 pollen record (35-18 cal ka BP) from Erhai Lake, southwestern China. The results disclosed a vegetation transition from temperate deciduous broadleaf forest dominance in late MIS3 to cold coniferous and mixed broadleaved/coniferous forests in LGM. This vegetation dynamic involved functional composition shifts from competitive-driven functional convergence to partial recovery via niche differentiation during the late MIS3, and finally to a low-diversity but functional differentiation state through trait complementarity and diversification strategies during the cold LGM. Our results likely support a function-mediated climate filtering process whereby climate change regulated long-term vegetation dynamics during the MIS3 to LGM transition primarily through shifts in CWM functional composition. These findings underscore the potential of pollen-based trait approaches to reconstruct ecosystem properties and advance our understanding of ecosystem change over decadal to millennial time-scales.

Climate change is expected to alter plant populations not only through direct environmental shifts but also via changes in biotic interactions, such as with herbivores and pathogens. As plant species are also expected to differ in their responses to both climate and antagonists, plant responses to both factors are expected to be variable and species-specific. To assess whether interactive effects of climate and antagonists on plant population dynamics are common and whether the strength and direction of plant responses vary across species, we conducted a multi-year field experiment that manipulated realistic climate change and experimentally reduced insect herbivores and fungal pathogens. We measured responses of plant vital rates, such as survivorship, growth, and reproduction across six grassland species. Using Integral Projection Models (IPMs) and Life Table Response Experiments (LTREs), we quantified changes in population growth rate across experimental treatments and the contribution of each vital rate to that observed change. Two of the study species declined so drastically over the course of the experiment that demographic quantification of population growth rates was not possible. From the remaining species, Bromus erectus and Plantago lanceolata show significant interactive responses of climate and antagonist reduction on population growth rates. In contrast, Dianthus carthusianorum and Tragopogon orientalis showed limited responses to experimental treatments. Notably, our results indicate that in some species biotic interactions may amplify the effects of climate change: the presence of plant antagonists exacerbates the negative effects of the future climate treatment on plant population dynamics. Our findings highlight the complexity in predicting plant population responses to climate change and provide insights for grassland management under future environmental conditions.